Draft version March 4, 2013 

Preprint typeset using LAT^X style cmulateapj v. 5/25/10 



OPTICAL AND INFRARED ANALYSIS OF TYPE II SN 2006BC 

Joseph S. Gallagher 1 , B.E.K. Sugerman 3 , Geoffrey C. Clayton 2 , J.E. Andrews 2 , J. Clem 2 , M.J. Barlow 4 , B. 
Ercolano 5 , J. Fabbri 4 , M. Otsuka 6-7 , R. Wesson 4 , and M. Meixner 7 
Draft version March i, 2013 

ABSTRACT 

We present nebular phase optical imaging and spectroscopy and near/mid-IR imaging of the Type 
II SN 2006bc. Observations reveal the central wavelength of the symmetric Ha line profile to be 
red-shifted with respect to the host galaxy Ha emission by day 325. Such an phenomenon has been 
argued to result from an asymmetric explosion in the iron-peak elements resulting in a larger mass of 
56 Ni and higher excitation of hydrogen on the far side of the SN explosion. We also observe a gradual 
blue-shifting of this Ha peak which is indicative of dust formation in the ejecta. Although showing a 
normal peak brightness, V ~ -17.2, for a core-collapse SN, 2006bc fades by ^6 mag during the first 
400 days suggesting either a relatively low 56 Ni yield, an increase in extinction due to new dust, or 
both. A short duration flattening of the light curve is observed from day 416 to day 541 suggesting 
an optical light echo. Based on the narrow time window of this echo, we discuss implications on the 
location and geometry of the reflecting ISM. With our radiative transfer models, we find an upper 
limit of 2 x 10 -3 M Q of dust around SN 2006bc. In the event that all of this dust were formed during 
the SN explosion, this quantity of dust is still several orders of magnitude lower than that needed to 
explain the large quantities of dust observed in the early universe. 

Subject headings: circumstellar matter — supernovae: individual (SN 2006bc) — dust, extinction 



1. INTRODUCTION 

The majority of dust formed in galaxies at the current 
epoch originates in the stellar wi nds of low m ass, asymp- 
totic giant branch (AGB) stars (|Gehrzlll989l) . Although 
the effective temperature of AGB stars, around 3000K, 
exceeds the condensation temperature for most composi- 
tions of dust grains, pulsation or convection driven shock 
waves propagating radially from the star can create tran- 
sitory dense, cool layers that are suitable for dust forma- 
tion. Once created the newly formed dust can be ushered 
away into the interstellar medium (ISM) via radiation 
pressure. 

Recent sub-mm observations of high redshift quasi- 
stellar objects (QSO's) and galaxi es have reveale d large 
quantities of du st at z > 5 (|Bertoldi et al.l 120031 : 
iDwek et"ahl 120071 ). Such redshifts correspond to an 
age of the Universe of 1 Gyr or less, considerably 
younger than the e xpected delay time for dust forma- 
tion in AGB stars (|Michalowski et alJfeOlOf ). An alter- 
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nate mechanism for the creation of this high redshift 
dust is dust formation in the ejecta of supernovae (SNe) 
given that the delay time associated with SNe is mea- 
sured in millions, not billions, of years . This was first 
sugges ted as a possible mechanism by iCernuschi et al.1 
(1967). Recent theoretical models predict that 0.1- 
1 M Q of dust should form in the ejecta of every SN 
event (|Todini fc Ferraral 120011: iNozawa et ail I2003D ; a 
production rate that can certainly account for the large 
amount of dust at high redshift (iMorgan fc Edmunds! 
[2001 iMaiolino et al.l 12001 IDwek et al.ll2007D . Observa- 
tional evidence of SN dust formation was first seen within 
the ej ecta of SN 1987A dLucv et al.l 119911: iMeikle et all 
119931 lErcolano et al.ll2007f) . This evidence came in the 
form of (1) the development of an infrared excess, (2) an 
increase in the optical extinction, and (3) the appearance 
of blue-shifted emission line profiles due to the greater at- 
tenuation of emission from the far-side, red-shifted ejecta 
emission. 

To date, evidence of dust formation has been seen in all 
types of core-collapse SNe, save for type Ic, and no signs 
of dust formation have been observed in thermonuclear, 
SNe la. In addition to SN 1987A, dust has been shown 
to ha ve condensed in the ej ecta of fi ve Type II SNe, 
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Kotak et al.1 l2009t lEabbri et al l 120111: ISueerman et al.1 
201H), and 2005af (jKotak et all 120061 ). In each case the 
quantity of dust formed was less than 10~ 3 M Q , con- 
siderably less than that needed to explain the dust ob- 
served at high-z. More recently, observations of SNe ex- 
hibiting strong circumstellar interaction have revealed 
a secondary location, outside of the ejecta, of dust 
formation. These observations have shown evidence 
of dust formation within the cool, dense shell (CDS, 
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iDeneault et al.l (|2003D ) between the forward and reverse 
shocks formed by the interaction of the ejecta and cir- 
cumstellar material (CSM). This mechanism h as been 
reported in four SNe. 1998S (iPozzo et al.1 l200l. 2005ip 
(iSmith et al.ll200l [FbTet al.ll2009fl. 2006ic (ISmith et al.1 
120081 iMattila et al.ll2008D. 2007it (Andrews et al.ll2011h . 
and 2007od ((Andrews et al.l 120101^ However, despite 
their potential to produce dust in both their ejecta and in 
their CDS, CSM-interacting SNe have been shown only 
to form comparable amounts of dust to their CSM defi- 
cient counterparts. 

Recent studies that have reevaluated the delay time 
limits for dust formation by AGB stars have concluded 
that AGB stars may be more efficient dust produc- 
ers at high redshift than first expected, but they have 
not been able to ru l e out SNe as major con tributers 
(jValiante et al.ll2009t iDwek fc Cherchnefil 1201 1ft . Other 
possible, non-stellar, sources of dust have been put for- 
ward s uch as the growth of dust grains in the ISM. lDraind 
( 2009) suggested that although SNe contribute a portion 
of the dust found in the ISM, that most of the dust mass 
observed in the ISM comes from the growth of these seed 
grains after they arrive in the ISM. 

In spite of the great interest in understanding the ori- 
gins of dust at high redshift, there is still no firm conclu- 
sion yet as to the ability of SNe to significantly contribute 
to the overall dust budget in the galaxy. The number of 
SNe that have been observed at multiple epochs and in 
multiple wavelengths is still very small, and the more SNe 
that are observed, the greater the realization that each 
SN is unique. This makes it of paramount importance 
to increase the sample of observed core-collapse SNe, so 
that a statistically significant conclusion can be made on 
their ability to form dust. To this end, we report multi- 
temporal, multi-wavelength observations of core-collapse 
SN 2006bc. In §2 we will describe our observing pro- 
gram and our reduction methods. In §3 we will report 
our results, and we summarize our conclusions in §4. 

2. OBSERVATIONS & REDUCTIONS 
2.1. Previous Observations and Relevant Statistics 

SN 2006bc was discovered on 2006 March 24.65 UT 
approximately 49" west and 43" north of th e center of 
its host, NGC 2397 (F igure m iMartinl d200l ). An opti- 
cal spectrum taken bv lPatat et al.l (|2006l ) with the ESO 
Very Large Telescope on 2006 March 28.06 UT showed 
narrow Balmer emission superimposed on a blue con- 
tinuum. The width of these Balmer lines suggested an 
expansion velocity of ^2000 km s _1 . Unresolved Na 
I D lines originating within the Milky Way and the 
host galaxy were also present, and both lines had mea- 
sured equiv al ent w idths of ~ 1.0 A (jPatat et al.1 120061 ). 
iPatat et"al~l ((2006) also determined a B— V color for 
SN 2006bc of ~ 0.2 while coincident Swift Ultravio- 
let/ Optical (UVOT) observ ations found U-B - -0.9 ± 
0.1 (llmmler fc Brownll2006ft . Based o n the color-SN type 
diagnostics of IPoznanski XalllpOOl . SN 2006bc is con- 
siderably bluer than most type II SNe with colors closer 
to that of type la SNe. However, the presence of hydro- 
gen emission and the similarities of its colors to type II-P 
SN 2 005cs and 2006at led to its classification as a type 
II-P (jlmmler fc Brownll2006D . 

SN 2006bc was further imaged on 2006 October 14 



(203 d) with the ACS Wide Field Camera (WFC) in the 
F435W , F555W, and F8UW filters by ISmartt et~aTl 
(2009). Alignment of these images with pre-explosion 
images taken with the Wide Field Planetary Camera 2 
(WFPC2) onboard the Hubble Space Telescope (HST) on 
2001 November 17 resulted in no detection of the progen- 
itor star.. However, based on the detection limit of the 
pre-explosion ima ges, and assuming th e progenitor to be 
a red supergiant, ISmartt et al.l (|2009l ) derived an upper 
mass limit for the progenitor of SN 2006bc of 12 M Q . 

Based on the equivalent width measurements of the 
unresolved Na I D lines for the NGC 2396 and the Milky 
Way by (IPatat et al.ll2006D. a nd using with the empir- 
ical relation of Bar bon et al.l (fl990l) relating EW(Na I 
D) with the intrinsic color excess, E(B-V), and R(V) = 
3.1, we find a substantial total foregr ound extinction t o 
SN 2006bc of Aw - 1.6. According to (jTullv et al.ll2"o"09h . 
NGC 2397 is at a distance of 20.3 Mpc, and we calculated 
the recessional velocity of NGC 2397 via redshifted HII 
emission lines in our own optical spectra to be (1462±12) 
km s~ \ Finally, based on negative detections bv lMartinl 
(2006) as late as 2006 March 16.59 UT, we adopt an 
explosion date that coincides with the date of discovery 
(JD 2453819.15). 

2.2. Optical Spectroscopy 

For this program data were taken with the Gemini 
South telescope utilizing the GMOS-S instrument dur- 
ing the 2007A and 2007B semesters (GS-2007A-Q-4 & 
GS-2007B-Q-39). Dates of observation include 2007 
February 12 (324 d), March 18 (358 d), May 11 (412 
d), September 20/24 (~544 d), November 11/14 (~596 
d), and 2008 January 16 (662 d; see Table 1). Dur- 
ing the first three epochs, 3 x 900s spectra were obtained 
in long slit mode. The exposure time was augmented 
to 3 x 1050s for the last three epochs to compensate for 
the expected SN optical decline. The GMOS-S spectro- 
graph, equipped with the B600-G5323 grating, allowed 
for 2760 A coverage and a FWHM resolution of ~3.9 
A. The central wavelengths of individual spectra were 
5950A, 5970A, and 5990A, respectively, to allow for ef- 
fective removal of chip gaps and to avoid having the gaps 
fall on important emission features. A 0.75" slit width 
was used with the B600-G5323 grating in first order along 
with 2x2 binning in the low gain setting. 

Spectra were reduced using the TRAJ^gemini package. 
An example of the slit orientation from day 324 is shown 
in Figure [2l The inset displays a section of the 2-D spec- 
trum including emission from Ha and [Nil]. The figure 
highlights the contamination of galactic emission in the 
light from the SN. The SN and sky regions used during 
the 1-D spectral extraction were determined by visual 
inspection. Care was taken when defining the regions of 
sky to include continuum emission from the galaxy, but 
to avoid stars and clumpy HII regions falling on the slit 
as shown in Figure [2l The 3 spectra obtained each indi- 
vidual night were averaged to improve signal-to-noise. 

2.3. Optical Imaging 

8 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Coincident optical imaging was obtained using the 
GMOS-S imager in three SDSS filters, cf, r 7 , and i. 
The images were reduced and stacked using the stan- 
dard routines within the gemini package. To remove as 
much contamination from galactic light as possible, we 
performed PSF-matched difference i maging in each filter 
using the IRAF difimphot package (|Tomanev fc Crottsl 
ll996tlSu"german eta l. 2005). The procedure entailed reg- 
istering and matching the PSFs of an object and refer- 
ence image. Given the lack of pre-explosion GMOS-S 
images of NGC 2397, we opted to use our day 662 </rY 
images as the reference epoch. In each filter, the refer- 
ence image was subtracted from the images of the earlier 
epochs, resulting in a final PSF-matched difference image 
in g'r'i' for days 324, 358, 412, 544, and 596. An example 
is shown in Figure [3] of our day 324 i' subtraction. 

PSF-photometry was then performed on the differ- 
ence images using the IRAF daophot package. The 
instrumental g'r'i' magnitudes were converted into the 
standard Johnson-Cousins V RI C using a transformation 
involving a least-s quares fit with a floating zero-point 
(jWelch et al.ll2007h . A BVRI C photometric sequence of 
secondary standards was derived for the 2006bc field. 
Observations of the se condary s tanda rd stars and ^30 
stars from the lists of lLandoltl ()2009D were taken with 
the Y4KCam CCD on the Yale 1.0- m telescope operat- 
ing at the Cerro Tololo Inter- American Observatory on 
2008 October 28. For full details on the data reduction 
and the cre ation of our photome tric sequence see Ap- 
pendix A of lAndrews et all (|2010| ). The standard stars 
are shown in Figure Q] and their corresponding magni- 
tudes are reported in Table 4 . PSF-photometry on the 
non-differenced images of these secondary standard stars 
allowed us to derive the unique, filter-dependent zero- 
points for each night. It is important to note that any 
remaining SN flux present in our reference epoch will 
lead to our underestimating the flux at a given epoch. 

Late time imaging data were obtained by HST using 
WFPC2 on 2007 September 09 (532 d) and 2008 Febru- 
ary 17 (694 d). Data were taken using the _F606VFand 
F814W filters, and the SN was centered on the Plane- 
tary Camera (PC) CCD chip. The images were pipeline 
reduced, cosmic rays were removed, and images were 
drizzled to a higher spatial resolution, 0.03" pixel -1 
(65% of the native pixel scale), using the STSDAS pack- 
age. Instrumental magnitudes were determined via PSF- 
photometry using the daophot package. We generated 
simulated PSFs for t his analysis using th e PSF creation 
software Tiny Tim (jKrist fc Hook|[200l . The F6061U 
and F814W instrumental magnitudes were transformed 
into Johnson V and I c , respectively, and charge trans- 
fer efficiency (CTE) correctio ns were applied using the 
iterative procedure detailed in iDolph in (2000l |2009( ) 

2.4. Near-Infrared Imaging 

Wide band F110W, F160W, and F205W images of 
2006bc were obtained on 2007 September 17 (541 d) 
and 2008 February 17 (694 d) using HST equipped with 
the Near Infrared Camera and Multi-Object Spectro- 
graph (NICMOS, Meixner: G011229). The NIC2 cam- 
era employing the square wave dither pattern was uti- 
lized. The images were pipeline reduced and photometry 
was carried out in similar fashion to that of the WFPC2 
data. Model PSFs were created with Tiny Tim and were 



used along with daophot to calculate instrumental mag- 
nitudes. HST system magnitudes (ST) were calculated 
from the measured count rate (CR) and the header key- 
words PHOTFLAM according to Equation (1). 

STMAG = -2.5 x log 10 {PHOTFLAM x CR) - 21.1 (1) 

PHOTFLAM is defined as the mean flux density in 
units of erg cm -2 s -1 A -1 that produces 1 count per 
second in a given HST observing mode. The ST mag- 
nitudes were then transformed into JHK magnitudes 
using the IRAF synphot package. 

2.5. Mid- Infrared Imaging 

Three epochs of Spitzer/IRAC (3.6, 4.5, 5.8, and 8.0 
fim) and MIPS (24/im) images were taken along with two 
epochs of IRS Peak-Up image (PUI) data at 16/im. The 
IRAC data were taken on 2007 September 12 (537 d), 
2008 February 12 (690 d), and 2008 May 08 (776 d). The 
IRS PUI data were taken on 2007 August 02 (496 d) and 
2008 January 01 (648 d) while the MIPS data were taken 
on 2007 September 15 (540 d), 2008 February 10 (688 d), 
and 2008 April 15 (753 d). The chronology of all observa- 
tions is summarized in Table 1. To improve photometric 
quality, all the data were mosai cked and resampled usin g 
the software package MOPEX (jMakovoz fc Kharil20l)5T ). 
The IRAC images were drizzled to 0.75" pixel" 1 (63% of 
the native pixel scale), the IRS PUI 16/xm to 0.9" pixel -1 
(50%), and the MIPS 24^111 to 1.2" pixel -1 (50%). Dif- 
ference imaging analysis was performed as described 
above with the day 776 (IRAC), day 753 (MIPS), and day 
688 (IRS /PUI) images subtracted from the earlier epochs 
for the three instruments, respectively. PSF-photometry 
was performed at the SN position using DAOPHOT 
while a customized implementa tion of the IRAF packag e 
DIMPHOT, called DIPHOT (|Tomanev fc Crottsl [19961) . 
was used to obtain photometric errors or upper detection 
limits for the differenced images. 

3. RESULTS 

3.1. Optical Emission Line Profiles 

As noted earlier, the spectroscopic analysis of SN 
2006bc is difficult, particularly at late times, due to its 
location within the bright SBb host galaxy, NGC 2397. 
Our six spectra, stretching from day 325 to 663, are 
shown in Figured The spectra exhibit prominent narrow 
emission lines of H/3, [OIII], [Nil], [SII], and Ha typically 
observed in the HII regions of star-forming galaxies. The 
spectra also show lines of [OI] and Ha with widths con- 
sistent with emission from the evolving SN ejecta (500- 
1000 km s -1 ). The evolution of [OI] and Ha are shown 
in Figures [5] and HI respectively. Although Figure [5] has 
been corrected for the velocity of the host galaxy, Figure 
[6] has not. Little change is observed in the [OI] emission 
profiles; however, in spite of strong contamination from 
nebular emission, noticeable evolution is observed in the 
profile of Ha. Figure [6] shows the evolution of the Ha + 
[Nil] emission profile for which we have carefully aligned 
the narrow host galaxy Ha and [Nil] lines. Apparent 
is the underlying emission from the SN fading over the 
course of our observing period. 
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In order to separate the nebular emission from the SN 
emission, we simultaneously fitted Gaussian curves to the 
blended Ha + [Nil] profile of the host galaxy and SN. 
The results of this analysis are given in Figure [3 The top 
two panels show the fit to our day 325 profile on the left 
and its Gaussian components on the right. Subsequent 
plots show the fits to our remaining five epochs. The 
most obvious characteristic of the SN emission is that 
by day 325 the central wavelength of the SN emission 
is red-shifted with respect to the galactic Ha emission 
by ^800-900 km s _1 . A similar phenomenon was ob- 
served i n the optical spectra of SN1987A, 1999em, and 
2004di (IPhillips fc William^[ll9l1:lElmhamdi et al.H2003t 
iChugai et al.l 120051 ). In the case of SN 1987A, the emis- 
sion peaks of several elements, including hydrogen, were 
observed to be redshifted from 600-1000 km s _1 during 
the period 150 to 6 00 days after e xplosion. It was subse- 
quently argued by IChug ai (199lJ) that this redshift was 
due to an asymmetry in the 56 Ni ejecta. Just like SN 
2006bc, SN 1987A did not show this asymmetry in the 
[OI] emission lines , but it did show strong asy mmetry 
in Fe II and Ni I (jMeikle et all 119891 ) . This led iChugail 
(|1991|) to suggest that the asymmetries involved only the 
iron-peak elements with the redshifting of the Ha line 
being a secondary effect due to an increase in the exci- 
tation of hydrogen in positions where the mass of 56 Ni 
is higher. Given that we are observing the emission line 
profile asymmetries in Ha, but not [OI], it is likely that 
we are seeing a similar effect in the optical spectra of SN 
2006bc. 

Figure [8] shows the evolution of our Gaussian model to 
the SN Ha emission profile with time. The inset shows 
the central wavelength evolution of these models which 
reveals an apparent blueshift of the Ha peak from day 
325 to day 663. Figure [3] shows a blueshift of approxi- 
mately 200 km s _1 over the course of our observing pro- 
gram. Although it is a possibility the such a blueshift 
could be a further consequence of the apparent asymmet- 
rical explosion, blueshifted emission lines are also often 
the result of the formation of new dust in the ejecta re- 
sulting in the emission from the back-side of the explosion 
experiencing greater attenuation than emission from the 
front-side. Consequently, the apparent blueshift of the 
peak in Figure [8] could be a sign of dust formation in the 
ejecta of SN 2006bc. 

3.2. Radiative Transfer Modeling 

From our PSF-matched photometry we were success- 
ful in obtaining detections in all four IRAC channels on 
day 537, but we did not detect the SN on day 690 in 
any channel. Furthermore, we did not obtain detections 
in the MIPS or IRS/PUI data at any epoch. We com- 
bined our IRAC detections on day 537, our HST results 
from day 532, and upper limits at 16 /im (day 496) and 
24 pm (day 688) into a single SED for 2006c approxi- 
mately 550 days past maximum light. The SED of SN 
2006bc is shown in Figure [9j The near and mid-IR flux 
measurements are shown in Table 2) 

Even with mere upper limits beyond 10 /im, the IRAC 
results show evidence for strong IR excesses attributable 
to warm dust. The dotted-line in Figure |H] represents our 
blackbody fit to the SED. The data are well fitted by the 
sum of a 4900 K blackbody, representing the emission 
from the hot, optically thick, ejecta, and a 500 K modi- 



fied blackbody due to the emission from a spherical dis- 
tribution of warm dust with an outer radius of ^570 AU. 
The modified blackbody is a normal blackbody subject 
to a A -1 emissivity law. 

Quantitative modeling of the dust around SN 2006bc 
proceede d using our 3D Mon te Carlo RT code MO- 
CASSIN ()Ercolano et al.l (|2005[) . and references therein). 
Although we have presented evidence that 2006bc ex- 
ploded asymmetrically, we do not have any strong obser- 
vational constraints on the degree of asymmetry. Con- 
sequently, the dust around SN 2006bc was modeled as 
a spherically, expanding shell with an inner radius, Ri n , 
and an outer radius, R ou t- We modeled the dust around 
SN 2006bc in two separate ways. The first we refer to 
as our "smooth" RT model. Since the source of the lu- 
minosity is the ejecta, our smooth model assumes that 
the source luminosity is spread throughout the shell and 
proportional to the local density. We further assumed a 
uniform dust distribution throughout the shell according 
to an r~ 2 density profile. For our second model, we allow 
the dust to form into clumps. For this reason, we call it 
our "clumpy" model. The clumpy model comprises an 
inhomogeneous distribution of spherical clumps embed- 
ded within an interclump medium of a supplied density. 
The irradiating photons are only produced in this in- 
terclump medium with the clumps assumed to be dark. 
The siz e of the spherical clumps are given by r r = 0.025 
(rout) (|Herant fc Wooslevl Il994t ISugerman etall I2006T ) 
with a volume filling factor, f c = 0.20, and density con- 
trast a = p c /p = 1x10 s . For a more detailed discussion 
of modeling a c l umpy environment with MOCASSIN see 
lErcolano et al.l (|2007f) . For each model, we adopted a 
standard grai n size distribution of a~ 3 5 between 0.005 
and 0.05 /im (ISugerman et al.l [20061: iMeikle et al.l [20071: 
iKotak et al.l I2009D . We tested two dust compositions 
with varying amounts of silicate (Si) and amorphous car- 
bon (AC) dust. The first was the AC dominated dust 
(75% AC, 25% Si), the second was Si dominated (25% 
AC, 75% Si ). Both models used th e optical gr ain con- 
stants from lOssenkopf et al.l (|1992l ) for Si and iHannerl 
(fl~98ll for AC. 

For both the smooth and clumpy models, the primary 
input parameters that we varied were the ejecta tem- 
perature and luminosity, the inner and outer radii of the 
dust shell, and the mass of dust present. Initial estimates 
of the ejecta luminosity and the extent of the dust shell, 
came from our optical photometry and the blackbody fits 
described above, respectively. These parameters, along 
with the dust mass, were then varied to produce the best 
fit to the optical/iJST and IRAC data points while si- 
multaneously trying to predict more flux than the upper 
limits at 16 and 24 /jm would allow. 

Our Si and AC dominated smooth models are both 
shown in Figure [9[ We found we were unable to predict 
sufficiently low flux at 24 pm for our silicate dominated 
models. This is not unexpected as silicate dust has strong 
emission features at 10 and 18 /im which acts to hold up 
the SED flux relative to that at 3-8 /zm. Our best fit 
model, therefore, corresponds to dust with a composi- 
tion of 75% AC and 25% silicate dust. The other best 
fit parameters are an ejecta temperature of 5500 K, an 
ejecta luminosity of 1.2 x 10 6 L , a dust shell inner and 
outer radius of 67 AU and 6700 AU, respectively, and a 



SN 2006bc 



5 



total dust mass of 4 x 10~ 4 M Q (see Table 3). This cor- 
responds to a A„ = 1.4. If we allow the dust to form in 
clumps, the total amount of dust becomes ~ 2.0 x 10~ 3 
Mq . Even if all of this dust formed during the supernova 
explosion, it would still be orders of magnitude less than 
a typical core-collpase SN would need to form in order to 
explain the amount of dust observed in the early universe 
(jBertoldi et al.|[2003HMorgan fc Edmunds! 120031 ) . 

3.3. Possible Light Echo 

The optical and near-IR lightcurve of SN 2006bc are 
presented in Figure [TOl After differenced imaging of our 
GMOS data, detections were made and magnitudes were 
measured on days 325, 359, and 413. The data on day 532 
and 694 are from HST/WFPC2. The photometry has 
not been corrected for foreground extinction. SN 2006bc 
undergoes a steep drop immediately following maximum 
light leading us to categorize 2006bc as a Type II-L SN. 
Assuming a total foreground extinction, Ay ~ 1.6 (see 
above), and a distance to NGC 2397 of 20.5 Mpc, SN 
2006bc reached a maximum brightness of V ~ -17.2. 
Although showing normal peak brightness for a core- 
collapse SN, 2006bc then fades by ^6 magnitudes during 
the first 400 days su ggesting a relatively low 56 Ni yield 
in the SN explosion (jKasen fe Wooslevll2009| ) and/or an 
increase in extinction due to new dust. Typically, such 
steep declines are seen in under- luminous SNe, though 
two Type II plateau SNe, 19 94W (|Sollerman et al.lll998ft 
and 2007od (j Andrews et al.|[2010h . were very highly lu- 
minous ( V < -18) during the photospheric stage and yet 
showed low late-time luminosities attributed to a combi- 
nation of low 56 Ni mass and dust extinction. 

By the onset of the nebular phase, SN 2006bc begins 
to decline with a rate in V and I c that is consistent with 
the radioactive decay rate 56 Co. However, we observe an 
approximate 100 day flattening of the lightcurve in both 
V and I c from days 416 to 541, after which a 56 Co de- 
cline rate resumes. The most likely explanation for this 
phenomenon is an optical light echo, and given the nar- 
row time window in which the flattening was observed, 
we can place constraints on both the distance and geom- 
etry of the circumstellar material (CSM) responsible for 
the echo. There are two obvious geometries that could 
explain the flattening of the lightcurve. The first one has 
the CSM lying in a torus surrounding the SN. In this 
case, due to light travel effects, the late onset and short 
duration of the flattening requires that the torus to be 
oriented close to face-on and have a thickness of no more 
than ~ 100 lt-day. Furthermore, the fact that flattening 
begins around day 413 suggests that the inner edge of the 
torus or ring must lie approximately 400 lt-day (1.04el8 
cm) from the geometric center of the explosion. This is 
consistent with the size of the main ring around 1987A. 

Another possible scenario is that the CSM is located 
within a thin spheroidal shell with a radius of approx- 
imately 400-500 lt-day. Given the time it would take 
for light to travel to all parts of the shell and then be 
scattered in the direction to Earth, one is likely to as- 
sume that such a geometry could not produce the short 
duration flatteni ng th at we have observed. However, 
Su german et alJ (|2005h have shown the following propor- 
tion to be an empirically close fit to the average scatter- 
ing function for MRJN0 dust over scattering angles be- 

9 IMathis et alJ (p"977f l 



tween 6 = to ~ 160°. The echo brightness scales as, 

E.B. oc [1.0 + 2.21 x 10~ 3 x (<9) 2 ] ~ 1-123 . (2) 

This relation indicated that the echo brightness falls off 
sharply with increasing scattering angle, 6. With a thin 
shell of radius 400-500 lt-day, the scattered light from 
the CSM at low 9 will begin to flatten out the lightcurve 
400 to 500 days after the explosion while the material 
at increasing 6 could continue to hold the lightcurve up 
for the subsequent 100 days. Either of these scenarios 
are possible, though the latter seems more likely as the 
former requires such a specific inclination (i > 80) of the 
torus to fit the data. 

3.4. Conclusions 

We have presented multi-temporal, multi-wavelength 
observations of the core-collapse SN 2006bc during its 
nebular phase. Based on an analysis of these observations 
we have drawn the following conclusions. 

1. An analysis of the evolution of the optical spec- 
tra of SN 2006bc shows evidence that the average 
central wavelength of the SN Ha emission profile is 
redshifted relative to the nebular emission from the 
galaxy between 325 to 663 days. This is likely due 
to the asymmetric nature of the explosion. The fact 
that this same redshift was not seen in the emis- 
sion lines of [OI] suggests the asymmetry was of 
the iron-peak elements only. Furthermore, we find 
evidence for a gradual blue-shifting of the Ha emis- 
sion line profile during our observing period due to 
the formation of new dust within the SN ejecta. 

2. Radiative transfer modeling of the dust around 
2006bc has revealed an upper limit of 4xl0" 4 M 
of new dust having formed during the SN explo- 
sion. If we allow the dust to form into clumps, our 
models predict as much as 2 x 10~ 3 M Q of new 
dust formed during the explosion of SN 2006bc. 
These results are consistent with the mounting ev- 
idence suggesting dust formation within the ejecta 
of core-collapse SNe is insufficient to explain the 
large quantities of dust observed in the early uni- 
verse. 

3. An approximate 100 day flattening of the light 
curve beginning ~ 400 days past explosion provides 
evidence for a light echo off a thin, spheroidal, shell 
of CSM surrounding SN 2006bc at a distance of 
400-500 lt-day. 
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Figure 1. Gemini/GMOS-S image of the field around NGC 2397 and SN 2006bc obtained 324 days after explosion with the i' filter. A 
higher resolution inset showing the immediate region around SN 2006bc is shown and was obtained with HST/ACS by Smartt ct al. (2009) 
as part of GO10498. The figure also shows the secondary standard stars that we used to transform from the instrumental to standard 
magnitude system (See Table 4). 
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Figure 2. GMOS-S i' image of SN 2006bc from 12 February 2007. The position of the SN is marked with a white arrow. The orientation 
of the slit for the February spectroscopic observations is depicted. The image shows the crowded environment in which SN 2006bc is found, 
and the inset shows the section of the GMOS-S spectrum around Ha. The respective origins of three emission regions along the spatial 
axis are marked with arrows in the image. Care was taken to avoid emission from extraneous stars or HII regions (such as those shown 
above and below the SN position) in the section of sky used in the background subtraction. 




Figure 3. Example PSF-matched differenced image from our CMOS data. Left: Day 324 i' image of SN 2006bc. Middle: Day 662 i' 
image of SN 2006bc. Right: Final subtracted image showing the SN with the nebular contamination considerably reduced. 
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Figure 4. Gemini/GMOS-S spectra SN 2006bc spanning 6 epochs. Data have not been corrected for redshift or reddening. 
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Figure 5. Evolution of the [OI] A6300,6363 emission lines across our six epochs. The spectra have been redshift corrected, but they have 
not been dereddend. 
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Figure 6. Evolution of the Ha + [Nil] emission line profile of SN 2006bc. Full profile shows strong contamination from narrow nebular 
emission. This emission, originating in the host galaxy, is labeled as "host." Broad-width emission originating in the SN ejecta is also 
labeled. 
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Figure 7. Gaussian fits to the respective Ha + [Nil] profiles observed at each epoch. The fits are shown in black while the data are shown 
in grey. The fits are comprised of four Gaussian curves, one for each the SN ejecta emission, the galactic [Nil] A6549 and A6584 emission, 
and the galactic Ha emission. As an example, the four constituent parts of the fit for day 325 are given in the top-right panel. The lower 
panels show the fits for subsequent epochs. 
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Figure 8. Evolution of the SN component to our Ha + [Nil] fits from Figure [7] Inset shows the central recessional velocity of the SN 
emission versus time. The plot reveals an apparent blucshift of ~200 km s — 1 of the SN emission profile from day 325 to 663. 
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Figure 9. Blackbody and MOCASSIN fits to the SED of SN 2006bc approximately 550 days past maximum light. The optical data were 
taken on day 532 with #ST/WFPC2. The JHK fluxes come from HST /NICMOS data taken on day 540. The 3.0, 4.5, 5.8, and 8.0 (im 
fluxes come from Spitzer /IRAC data taken on day 537, while the 16 and 24 (im points were obtained with Spitzer using IRS/PUI (day 
496) and MIPS (day 688), respectively. The dot-dash line is best blackbody fit to the data, the solid line is our best fit AC dominated, 
"smooth" MOCASSIN radiative transfer model, and the dashed line is our best fit Si dominated, "smooth" model of the dust around SN 
2006bc. The data reveals a strong infrared excess by day 550. Our AC dominated model represents the better fit to the data (particularly 
at 24 /i), and predicts 4 x 10 -4 Mq of dust around SN 2006bc. If we allow for the dust to form into clumps, the total quantity of dust 
increases to 2 x 10~ 3 Mq. 
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Figure 10. VIcJHK lightcurve of SN 2006bc. The data within the first 10 days are unfiltered discovery magnitudes while the data 
be tween day 40 and 60 comes from Astrosurf.com. The V and I data points from day 203 are measured from the ACS WFC images taken 
bv lSmartt et all 120091 ). The data beyond 300 days comes from the observations made for our study. The lightcurves show an abrupt drop 
of ~7-8 magnitude during the first 400 days past maximum light. This is likely due to a combination of a low 56 Ni yield in the explosion 
and an increased optical extinction due to dust formation. The data also show a short-lived flattening in the light curve from ~400 to 500 
days past maximum due to an optical light echo from CSM. 
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Table 1 

Summary of Observations 



Epoch a 


JD 


Telescope 


Instrument 


Filters 


Exp. Time (s) 


203 


2454023 


HST 


ACS/WFC 


F555W; F814W 


1500; 1600 


3 21 


O A K A "1 A A 

2454144 


Gemini South 


GMOS-S spectra 




3x900 


324 


2454144 


Gemini South 


GMOb-S imaging 




60 


358 


O A K A "1 TO 

2454178 


Gemini South 


GMOS-S spectra 




3x900 


358 


2454178 


Gemini South 


GMOS-S imaging 


/ _/ 7 


60 


412 


O A F ,1000 

2454232 


Gemini South 


GMOS-S spectra 




3x900 


412 


O A F ,1000 

2454232 


Gemini South 


GMOS-S imaging 




60 


496 


2454olo 


Spitzer 


TTDC /DTTT 

IHb/FUl 


Blue 


31 


532 


O ^ F A O F O 

2454353 


HST 


Wt PC2 


b 606W; b 814W 


1600; 1600 


537 


2454356 


Spitzer 


IRAC 


CH 1-4 


10 


540 


2454359 


Spitzer 


MIPS 


24 fivcL 


31 


541 


2454361 


HST 


JN1CMOS 


tllOW; tl60W; t205W 


640; 512; 576 


544 


2454364 


Gemini South 


GMOS-S spectra 




3x1050 


544 


2454364 


Gemini South 


GMOS-S imaging 


gVi' 


60 


596 


2454416 


Gemini South 


GMOS-S spectra 




3x1050 


596 


2454416 


Gemini South 


GMOS-S imaging 




60 


662 


2454482 


Gemini South 


GMOS-S spectra 




3x1050 


662 


2454482 


Gemini South 


GMOS-S imaging 


gVi' 


60 


648 


2454467 


Spitzer 


IRS/PUI 


Blue 


31 


688 


2454507 


Spitzer 


MIPS 


24 /im 


31 


690 


2454509 


Spitzer 


IRAC 


CH 1-4 


10 


694 


2454514 


HST 


WFPC2 


F606W; F814W 


1600; 1600 


753 


2454572 


Spitzer 


MIPS 




31 


776 


2454595 


Spitzer 


IRAC 


CH 1-4 


27 



a Assuming explosion date JD = 2 453819 

b Data taken as part of GO10498 ()Smartt et al. 2009). Data were acquired from the archive, re-reduced and analyzed. 



Table 2 

Near and Mid-IR fluxes of SN 2006bc for day ~ 550 

A (/J,m) Flux (pJy) 

1.2 4.9+g;J 

1.6 5-7t?;o 

2.2 7.7±o7 

3.6 17.3±lg;l 

4.5 19.0±lg;| 

5.8 re.itHi 

8.0 145.0^239 

16.0 < 163.0 ' 

24.0 < 62.0 
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Table 3 

Best Fit Monte Carlo Radiative Transfer Models 



Model 


AC/Si 


T (K) 


Rin (AU) 


R oui (AU) 


L (L Q ) 


Tv 


M d (M Q ) 


Smooth 


0.75/0.25 


5500 


67 


6700 


1.2c6 


1.27 


4e-4 


Clumpy 


0.75/0.25 


5500 


67 


2000 


1.2c6 


1.00 


2c-3 



Table 4 

Secondary Standard Stars 



ID V (<5V) R (<5R) I (51) 



A 17.140 (0.009) 16.541 (0.006) 16.032 (0.023) 

B 17.003 (0.008) 16.367 (0.010) 15.826 (0.011) 

C 16.548 (0.013) 16.118 (0.009) 15.723 (0.051) 

D 16.878 (0.012) 16.180 (0.013) 15.624 (0.017) 

E 17.618 (0.011) 17.050 (0.009) 16.539 (0.106) 

F 16.348 (0.011) 16.014 (0.005) 15.656 (0.010) 



